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ABSTRACT 

We use sequences of images and magnetograms from Hinode to study mag- 
netic elements in internetwork parts of the quiet solar photosphere. Visual in- 
spection shows the existence of many long-lived (several hours) structures that 
interact frequently, and may migrate over distances ~ 7 Mm over a period of 
a few hours. About a fifth of the elements have an associated bright point in 
G-band or CallH intensity. We apply a hysteresis-based algorithm to identify 
elements. The algorithm is able to track elements for about 10 min on average. 
Elements intermittently drop below the detection limit, though the associated 
flux apparently persists and often reappears some time later. We infer proper 
motions of elements from their successive positions, and find that they obey a 
Gaussian distribution with an rms of 1.57 ± 0.08 km/s. The apparent flows 
indicate a bias of about 0.2 km/s toward the network boundary. Elements of 
negative polarity show a higher bias than elements of positive polarity, perhaps 
as a result of to the dominant positive polarity of the network in the field of 
view, or because of increased mobility due to their smaller size. A preference for 
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motions in X is likely explained by higher supergranular flow in that direction. 
We search for emerging bipoles by grouping elements of opposite polarity that 
appear close together in space and time. We find no evidence supporting Joy's 
law at arcsecond scales. 

Subject headings: Sun: photosphere - Sun: granulation - Sun: magnetic fields 



1. Introduction 



Magnetic elements have been extensively studied in network that partiall y outlines the 

boun daries of supergran ular cells. They were first observed as "magnetic knots" (IBeckers fc Schroter 
19681 ) and as "fil igree" (IDunn fc Zirkerlll973f ) , befo re be ing reso l ved in to strings of adjacent 
bright points bv iMehltretterl Jl974T ). iMullerl (Il977l ) and Iwilsonl Jl98lh showed that faculae, 
filigree, and bright points i n wide-band CallH filtergrams are manifestations of the same 
phenomenon. IMullerl (119831 ) introduced the name "network bright point" , and subsequently 
initia ted extensive studies of magnetic elements as G-band bright points (IMuller fc Roudier 



1984). Studies of bright points u s ing h i gh-resolution imaging (e.g., IBerger et al 



2007a 



1998a lbl. 12004: IBerger & Title 



Ishikawa et al 



2007 



1996 



2001 



Berger et al 



Rouppe van der Voort et al. 12005: Rezaei et al. 



1995 



20071 ) have since established that network bright 



points are manifestations of small, kilogauss magnetic elements that form the magnetic net- 
work 



1972 



Chapman & Sheelev 


1968; 


Livingston & Harvev 


1969; 


Howard k Stenflolll972: 


Frazier & Stenflo 



Stenflolll973f ). 



Magnetic field in internetwork ha ye been largely ignored until recently. However, it is 



currently being studied vigorous 



2003 




Truiil 


2006 




2007 





Socas-Navarro et al. 2004 



Y (e.g. 



Dommguez Cerdena et al.l 20031; Sanchez Almeida et al 



Lites k Socas-Navarrol booi ISocas-Navarro fc Litesl I2004J: 



o Bueno et a 



Rezaei et a 



Litcs et al. 



.1120041: Manso Sainz et a 



2007bl; Sanchez Almeida 



2004 



2007 



Khomenko et al.ll2005 



Harvey et al 



2007 



Dommguez Cerdena et al 



Orozco Suarez et al. 



2008 



\) . Many of these studies focus on determining the strength and dis 



tribution of flux. While there is some disagreement between results, it seems that field is 
ubiquitously present in internetwork at small scales. 

Concentrations of flux that are suf ficiently s trong may form internetwork bright points. 
Their existence was already noted by Muller (|l983). Few recent studies have analyzed 



these internetwork magnetic eleme nts. ISanchez A 



bright point density and lifetime, iDe Wijn et al 



meida et al.l (120041 ) measured internetwork 



J2005h reported t hat internetwork bright 



points trace locations of flux that may persist for periods of hours, iTritschler et al.l (120071 ) 
analyzed morphology, dynamics, and evolution of bright points in Ca II K in quiet sun, and 
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Sanchez Almeida et al.l (120071 ) searched for photo-spheric foot points of transition region loops 
in quiet sun. 

Magnetic elements were first modeled a s "flux tubes" by ISpmitl ffl973). Over the 



years, models grew increasingly complex (e.g., Knol 



Solanki fc Brigljevidll992l ; iGrossmann-Doerth et al 



ter fc Schuss. 



1994 



1998 



er 



1988 



Kel 



er et al. 



Steiner et al 



1998 



1990 



ISteiner 



20051 ). The explanation of photosph eric brightness enhancement of facul ae due to ho t walls 
pro posed bvlSpruit fc Zwaanl (119811 ) was verified by MHD simulations by iKeller et al.l (120041 ) 
and lCarlsson et al.l (120041 ) . On disk, bright points are formed as a result of radiation escaping 
from deeper, hotter layers due to the fluxtube Wilson depression. 



line cells on mesogranular scales (e.g., 


Domfnguez Cerdena et al. 


2003; 


De Wiin et al. 


2005 


Tritschler et al. 


2007; 


Lites et al. 


2008 


), while 


Berger et al. 


dl998 


b) observed "voids" in ac- 



tive network. In additi on, recent simulations in dicate that field concentrates on boundaries 



of mesogranular cells (jStein fc Nordlundl 120061 ). One would expect such a pattern to be 



set by granular motion s, similar to supe rgranular flows that eventually advect internetwork 



field into network (e.g.. lLisle et al.ll2000l ). Perhaps magnetic eleme nts form these patterns a s 



a result of flows associated with "trees of fragm enting granule s " (jRoudier fc Mullen 12004 ) 



which were previously linked to mesogranules by iRoudier et al.l (120031 ) . Flux is expunged by 
the sideways expansion of granular cells, and is collected in the downflows in intergranular 
lanes. In a "tree of fragmenting granules" , these flows would be expected to drive flux not 
only to the edges of individual granules, but also to the edges of the tree. 

In this paper, we present a study of the dynamics of magnetic elements in internetwork 
parts of the solar photosphere. This study is motivated by its relevance to the operation of a 
turbulent granular dynamo, the nature of quiet-sun magnetism, the generation of MHD waves 
that may propagate into the transition region and corona, and the coupling of internetwork 
field to the magnetic network. First, examples of magnetic elements are discussed in the 
context of fluxtube dynamics and lifetime (Sect. [37T]) . Magnetic elements are compared with 
bright points in G-band and CallH intensity in Sect. 13.21 A feature-tracking algorithm is 
applied in order to analyze the lifetime (Sect. I3.3P and the dynamics (Sects. 13.41 and 13751) of 
magnetic elements. Finally, a search for emerging bipoles is presented in Sect. 13.61 



2. Observations and Data Reduction 



We use an image sequence o f a quiet area recorded by Hinode (IKosugi et al.l 120071) us- 
ing the Solar Optical Telescope (ITsuneta et al.l 120081 ; ISuematsu et al.l 120081 ; llchimoto et al. 
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Fig. 1. — A set of sample images from the sequences. Clockwise from top left: Fel intensity, 
G-band intensity, CallH intensity, and Fel magnetogram scaled between —500 (black) and 
+500 Mx/cm 2 (white). Solid white contours drawn on the magnetogram outline the network. 
The area contained by the dotted white lines is not affected by incomplete sampling due to 
image motion. Regions of interest examined in Figs. [21 [31 and |5] are outlined by black lines. 
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20081 : IShimizu et al.ll2007l : Matsuzaki et al.ll2007t ) from 00:18 to 06:00 UT on March 30, 2007 



Hinode was programmed to observe quiet sun near a small area of weak plage. The center 
of the field of view was at (490", 50") at the beginning of the sequence, and approximately 
followed solar rotation during the sequence. Because \x ranges from 0.91 to 0.75 over the 
field of view, care must be taken to correct measurements of position and velocity for fore- 
shortening. Image sequences were recorded in the G band and in the CallH line using the 
Broadband Filter Imager. The latter have chromospheric contributions, but mostly sample 
the upper photosphere due to the broad filter bandwidth of 0.3 nm. The Narrowband Filter 
Imager was used to record Stokes I and V in the photospheric Fel line at 630.2 nm at an 
offset of —12 pm from line center. The spectral resolution of the filter is 9 pm at this wave- 
length. The Stokes V signal is sensitive to Dopplershift as a result of flows along the line 
of sight. It is impossible to compensate for this effect due to the single line position used 
in the observations. The frames consist of 1024 x 512 square pixels. Each pixel in a frame 
corresponds to 2 x 2 pixels on the CCD, summed to sacrifice resolution for increased cadence 
and field of view at constant telemetry. The G-band and Ca II H images have a pixel scale of 
0.11" and a field of view of 112" x 56". The Fel images have a pixel scale of 0.16", resulting 
in a larger field of view of 164" x 82". A total of 582 frames were recorded in each passband 
at regular cadence of 35.14 s. 

The G-band, Ca II H and Fe I filtergrams were corrected for dark current and were flat- 
fielded using the SolarSoft procedure fg_prep. For the Fel data, a custom flatfield was 
derived from Fel images recorded between 00:18 and 09:25 UT on March 30, 2007. This 
comprises the sequence described above, and a slightly shorter, but otherwise identical se- 
quence of 324 frames of the same area. A synoptic observation separates the two sequences. 
The Fel flatfield is created by averaging the raw data of both sequences. Visual inspec- 
tion of the resulting flatfield shows little remaining signal of solar origin, and clearly shows 
the CCD fringe pattern. We apply a simp le calibration of the Stokes V images to produce 



magnetograms following IChae et al.l (120071 ) . 



We carefully align the frames using Fourier cross-correlation techniques. The sequences 
of magnetograms, G-band images, and Ca II H images are first aligned separately. The dis- 
placements computed from magnetograms are also applied to associated Fel filtergrams. 
Next, the G-band and CallH sequences are aligned to the Fel intensity and unsigned mag- 
netogram sequences, respectively. Figure [1] shows a sample Fe I filtergram, the associated 
magnetogram, G-band image, and CallH image. 

A mask of network areas is created by taking a suitable threshold of the average mag- 
netogram, smoothed to remove structures with sizes below 2". The network mask is shown 
in the bottom left panel of Fig. [IJ 
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Fig. 2. — A sample set of X—t slices of the magnetograms at successive Y positions. The 
region is indicated in Fig. [TJ Each X—t slice shows the density of flux along the line of sight 
averaged over 3 pixels (0.48") in Y at the indicated Y position. The step size in Y is also 3 
pixels. The flux density is scaled between —200 (black) and +200 Mx/cm 2 (white). 



We identify areas of significant magnetic flux along the line of sight using a hysteresis- 
based algorithm. In the remainder of this paper, "flux" refers to flux along the line of 
sight, unless otherwise indicated. Such an algorithm searches for features in data using a 
low threshold, but a ccepts a feature only if it is also detected using a high threshold (e.g., 



DeForest et al.ll2007l ). The low threshold must be chosen such that features close to the noise 
level are still identified, while the high threshold minimizes the number of false positives. 
We begin by producing a running average of three frames of the sequence of magnetograms. 
Then, each frame is convolved with a round kernel with a diameter of 5 pixels, and also 
with one with a 11-pixel diameter. At this point, we split the analysis between positive 
and negative flux by the sign of the frame convolved with the 11-pixel kernel. We apply 
a suitable threshold to the result of the 5-pixel convolution for the high-level filter. The 
low-level filter is constructed by application of a threshold on the difference of convolutions 
with the small and large kernel. Detections in the areas affected by incomplete sampling due 
to image motions are discarded in both filters. Several artifacts caused by hits of cosmic rays 
are also removed. Finally, a detection of the low threshold is accepted if there are detections 
reaching the high threshold at three or more different times. 

The algorithm identifies 11579 magnetic elements with positive flux, and 5226 with neg- 
ative flux. Visual inspection of the result shows that this procedure rejects noise adequately 
and often captures weak fields. Obviously, structures with very short lifetimes, i.e., less than 
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Fig. 3. — A second set of sample X—t slices of the magnetograms at successive Y positions. 
The format is identical to the format in Fig. [2j 



three time steps, are discarded by our algorithm. We believe this is no serious issue for the 
present analysis, because we are primarily interested in longer-lived structures. 

In this paper, we focus on small-scale concentrations of vertical field in internetwork 
areas, that we will refer to as "internetwork magnetic elements" or IMEs. We discard 
elements when they overlap with network areas. 



3. Analysis, Results, and Discussion 

3.1. Visual inspection 

We searched for patches of IMEs in the sequence of magnetograms. To this end, we 
employed a "cub e sheer" dissec t ing t he data cube into X-Y and X-t slices, similar to 



the procedure by iDe Wijn et al.l (120051 ). who inferred the existence of long-lived "magnetic 
patches" in quiet sun internetwork areas from recurrent bright points in G-band and Ca II H 
sequences. We found that while many long-lived structures of one or more IMEs exist, they 
do not typically show strong location memory on timescales of several hours. Some remain 
stationary for long periods, while others may migrate significant distances in a few hours. 
MEs often interact with other MEs while migrating, making the identification of well-defined 
structures difficult. In addition, many short-lived concentrations of flux appear frequently in 
the vicinity of longer-lived structures, but with unclear association to a particular structure. 
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We cannot combine IMEs into "patches" following the procedure of iDe Wijn et al. 



( 120051 ). Their proxy- magnetometry likely misses many weak and short-lived structures. In 
the present data, the density of IMEs, especially of short-lived flux concentrations, and the 
many interactions of long-lived structures inhibits successful grouping of elements into well- 
defined patches. The parameters that govern grouping of elements can be tuned to yield 
either many short-lived patches consisting of a single IME, or a small number of long-lived 
patches containing the vast majority of IMEs. 



Figure [2] imitates "cube slicing" (cf. Sect. 2 of lDe Wijn et al.ll2005l ) on one of the regions 
of interest indicated in Fig. [IJ It shows a sequence of X-t cutouts at progressive Y loca- 
tions. The left edge of the region is in the magnetic network (see Fig. [T|). Many short-lived 
concentrations of flux with both positive and negative polarity are visible throughout the 
region. In addition, several elements with negative polarity show up as trails toward net- 
work in the slices from Y = 30.9" to Y = 33.3". They start around X pa 465" at t = 0. At 
t ~ 150 min, their motion in X slows as they reach the stationary element around X pa 455" 
and Y pa 30.0". True cube slicing shows that while the IME is continuously present, it 
sometimes drops below the detection limit and reappears some time later. 

Figure [3] shows a different region in the same format as Fig. [2j In this case, an IME 
of positive polarity appears around X pa 460" at t ~ 90 min, and remains visible until the 
end of the sequence. It remains stationary at first, then migrates quickly in X starting from 
t pa 270 min. Again, many short-lived concentrations of flux of both polarities are visible 
throughout the region. Some short-lived elements appear to be part of larger structures. In 
the Y = 56.7" panel, for instance, the concentration of negative polarity around X pa 453" 
and t pa 270 min seems connected to the stable element around X pa 455" and t ~ 230 min 
by a faint trail of negative flux best visible in the Y = 56.2" panel. Similarly, intermittent 
concentrations of positive flux around X pa 454" and t pa 30 min at Y — 55.2" may be 
connected to the long-lived structure around X pa 459" and Y = 55.7", as may elements 
around X pa 455" and t pa 240 min at Y = 54.8". 

Figure H] shows a sample mask of magnetic elements. It corresponds to the frames 
shown in Fig. [TJ IMEs are not spread homogeneously over internetwork, but rather appear 
to outline cells on scales of several Mm. 



3.2. Comparison with bright points 



We search for bright points in the G-band and Ca II H image sequences at those locations 
where our algorithm detects an IME. A random IME at a random time during the sequence 
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Fig. 4. — Sample mask of magnetic elements corresponding to the frames shown in Fig. [TJ 
Elements with positive flux are blue, those with negative flux red. Network is shaded gray. 



is selected, and the magnetogram at that location is manually compared to the G-band and 
CallH intensities. Fragmented IMEs are treated as if they were multiple IMEs. IMEs that 
have no clear visually identifiable flux are discarded. Such IMEs typically have strong Stokes 
V signal in the preceding or following frame. They are detected because the identification 
algorithm averages over three time steps. We so inspect 400 IMEs. Figure [5] shows three 
samples. A bright point is identified in the G-band intensity in 88 cases (22%), and in 
the CallH intensity in 73 cases (18%). 67 cases (17%) exhibit a bright point in both the 
G-band and CallH intensities. Due to the small number of samples, the uncertainty on 
these measurements amounts to 5%. While the difference between G band and CallH is 
statistically insignificant, the numbers do agree with the impression that bright points are 
more easily visible in the G-band than in the Ca II H line. Weak bright points are drowned 
in the background of reversed granulation in CallH, while G-band bright points appear in 
the dark intergranular lanes. 



These findings agree with ground-based results of lBerger fc Title! (120011 ) and llshikawa et al. 



( 120071 ) . who found that magnetic field is a necessary, but not a sufficient condition for the 
formation of bright points in G-band or CallK intensity. The visual inspection reveals no 
evident relationship between morphology of the IME and the existence of a bright point in 
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the G-band or CallH intensity. Some IMEs appear small and weak, either because the field 
is weak, or because it is oriented away from the line of sight, yet have clear bright points 
(top row of Fig. [S]), whereas other strong elements do not show any brightening (bottom row 
of Fig. E]) . 

Most IMEs that have an associated Ca II H bright point also have a bright point in G- 
band intensity. This suggests that the mechanism of brightness en hancement is similar. G 



band bright points are formed by weakening of molecular CH lines (lUitenbroek fc Tritschlei 



2006) , as a result of parti al evacuation of the magnetic element ( jSpruit 



1976 



Keller et al. 



2004 ; ICarlsson et al.ll2004j ). In the wings of CallH & K, the magnetic element is c ooler at 
equal geometrical height, but hotter at equ al optical depth (ISheminova et al.l 120051 ) . Since 
the CallH filter of Hinode is 0.3 nm wide (ITsuneta et al.ll2008l ). the bulk of the emission 
is formed in the wings of the line. The root cause of brightness excess in both the G-band 
and CallH filtergrams used here is the fluxtube Wilson depression. The chromospheric 
contribution of the core of the CallH line will become more important as narrower filters 
are employed, and it is likely that correspondence between bright points in the G-band and 
Ca II H intensities will be reduced. 



3.3. Lifetime 



We next compute lifetimes for each IME. Figure [6] shows the histograms of IME life- 
times, for positive and negative polarity separately. Elements have an average measured 
lifetime of about 10 min. It should be n oted that accurately determining the lifetime is 
difficult and prone to errors (see Sect. 3 of lDeForest et al.l 120071 ). 



The average lifetime resembles results of IBerger et al.l (ll998bl ). who measured a mean 
lifetime of 9.3 min for bright points in network areas. However, it is significantly longer 
than those m easured for internetw ork bright points in G-band (3.5 min) and CallH images 
(4.3 min) by lDe Wijn et al.l (120051 ) . Based on positions of their bright points, those authors 
argued that bright points map positions of long-lived fields that exist before a bright point 
appears and persist after it disappears. 

The present analysis allows us to track magnetic elements that no longe r have an asso- 



ciated bright point, yet we do not find lifetimes of several hours, predicted by lDe Wijn et al. 



( 120051 ) based on statistical analysis of groups of bright points in a 1-hour sequence. However, 
it is clear that many elements live longer than our detection algorithm is able to track them 
from visual inspection of the magnetograms (cf. Sect. 13. ip . The lifetime computed here is a 
measure of how long flux remains concentrated enough for our algorithm to track it, rather 



Fe I magnetogram 



G-band intensity 



Ca II H intensity 




Fig. 5. — Sample IMEs with (top two rows) and without (bottom row) associated bright 
points. The IME mask is indicated by the contour in the magnetogram. The magnetogram 
is scaled between —200 (black) and +200 Mx/cm 2 (white). The top row shows an IME that 
has little apparent flux, yet still produces a clear bright point in both the G-band and Ca II H 
intensity. The IME in the bottom row is similar to the IME in the middle row, but does not 
have associated bright points in either the G-band or Ca II H intensity. 
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Fig. 6. — Histograms of the lifetime of IMEs. Blue: positive polarity. Red: negative polarity. 
The averages are indicated by the vertical dashed lines. 

than an estimate of the lifetime of flux itself. 



We measure the positions of IMEs by computing their center of mass, using the flux 
density from the magnetograms for weighing. First, elements are discarded if they reach 
the edges of the area affected by incomplete sampling due to image motions (dotted white 
lines in the bottom left Fig. [1]). Splitting or merging elements are then cut up into separate 
detections by wiping out the IME in the frame after the merger, or in the frame before the 
split. Finally, we identify elements whose mask changes between consecutive frames by over 
five times the minimum mask size in those frames. This can be caused by the appearance of 
a new element close to an existing element, so that it is identified as the existing element. 
For those cases, the IME in the second frame is erased. For the remaining detections, the 
center of mass of the mask is computed, while weighing with the magnetograms smoothed 
over three frames in time. We so find 136 271 and 76211 positions for IMEs with positive 
and negative flux, respectively. 

Hereafter we use the term 'velocity' for the inferred proper motions of IMEs. A velocity 
can only be computed if a magnetic element has a well-defined position in two consecutive 
frames. This reduces the number of usable positions to 112 871 measurements for IMEs with 
positive flux and 62 748 for those with negative flux. The top two panels in Fig. [7] show the 
distribution of proper motions in X and Y. They are nearly Gaussian around the origin with 
a slight overdensity at velocities close to zero and in the far wings. Disregarding polarity and 
direction, a Gaussian fit gives an average of 0.02 ±0.10 km/s and an rms of 1.57 ±0.08 km/s. 

The rms of the proper motions is nearly a factor of three higher than that measured 



3.4. Proper motions 
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Fig. 7. — Histograms of the velocity of IMEs in X (top left panel) and Y (top right panel), 
and of the speed of IMEs (bottom left panel). Blue: positive polarity. Red: negative 
polarity. The bottom right panel shows velocity in X vs velocity in Y, combining both 
polarities. Dashed lines indicate the first moments. 



by Ivan Ballegooijen et al.l (119981 ). who found an rms of 0.54 km/s using a cork tracking 
technique. However, t he present analy s is is i n good agreement with the measurement of an 
rms of 1.31 km/s from iNisenson et al.l (120031 ) . Those authors attribute the higher velocities 
to the selection of dynamic isolated bright points, rather than stable network elements. Since 
IMEs are typically also isolated, the results presented here are affected by the same selection 
effect. 

The histogram of the speed shows the expected Rayleigh distribution, though it has a 
noticeably longer tail toward higher speeds. Finally, there is no evidence of a correlation 
between the velocity in X and the velocity in Y in the scatter plot (Fig. lower right panel). 

In total, there are 2426 measurements (1.4%) of proper motions larger than 8 km/s. 
Investigation shows that these high velocities are caused by inaccuracies in the determination 
of the center of mass of weak elements as a result of noise in the magnetograms. 

Figure [S] shows the autocorrelation of the velocities in X and Y. Errors in measur ement 
of position propagate into the autocorrelation of the velocity (cf. INisenson et al.ll2003l ). The 
first four measurements around At = are influenced by the present analysis because it 
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Fig. 8. — Autocorrelations of the velocity of IMEs in X (left panel) and Y (right panel). 
Blue: positive polarity. Red: negative polarity. 



averages three frames in time. While this precludes determining the correlation time of IME 
velocities, the autocorrelations are consist ently positive fo r delay s up to 600 s, in agreement 
with results from bright point tracking by lNisenson et al.l (120031 ). 



3.5. Direction 

We speculate that the low, yet positive, autocorrelation at long delays is caused by a 
consistent, slow drift of the IMEs toward network. The component of velocity toward the 
nearest network boundary was computed over one time step, corresponding to 35 s. Elements 
that are closer to the area affected by image motion than to the nearest network pixel are 
discarded. The results are shown in Fig. [9j IMEs display a broad velocity distribution, 
with a slight bias toward the network boundary. Using a Gaussian fit, we measure a bias 
toward the network boundary of 0.152 ±0.006 and 0.218 ±0.009 km/s, and rms velocities of 
1.746 ±0.005 and 1.840 ±0.007 km/s, for IMEs with positive and negative flux, respectively. 

Over 17 time steps (597 s), close to a granular turn-over timescale, the average velocity 
changes only little, to become 0.104 ± 0.006 and 0.227 ± 0.010 km/s for IMEs with positive 
and negative flux, respectively. The rms is reduced significantly to respectively 0.608 ±0.005 
and 0.703 ± 0.008 km/s. 

The average velocity toward the network boundary is higher for IMEs with negative flux. 
The difference between the polarities is statistically significant up to velocity measurement 
over 35 time steps (1230 s), and tends to grow as the number of time steps increases. The 
number of measurements drops below acceptable levels for longer delays. Perhaps the higher 
mobility of elements with negative flux is a result of an attraction to the network, which 
consists predominantly of positive flux. Alternatively, it may be a result of higher mobility 
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Fig. 9. — Histograms of the velocity of IMEs toward the network boundary measured over 
one time step (35 s; left panel) and 17 time steps (597 s; right panel). Blue: positive 
polarity. Red: negative polarity. The best-fit bias velocities are indicated by the vertical 
dashed lines. Both polarities show a broad distribution with a slight bias toward velocities 
in the direction of the nearest network pixel. The width of the distribution diminishes with 
increasing timescales, while the bias does not change noticeably. 



of IMEs of negative polarity due to smaller size (median size of 14 pixels in the detection 
mask) compared to IMEs of positive polarity (median size of 17 pixels). 

The left panel of Fig. [10] shows the distribution of the directions of velocity over one 
time step. It shows a preference for velocities with directions around 0° and 180°, i.e., along 
the X axis. Correspondingly, the peak of the velocity distribution is higher in Y than it is in 
X (top panels in Fig. [7j), and the scatterplot of Vy against vx is somewhat oval, indicating 
there are more IMEs moving at higher speed in the X. In the classical view that IMEs are 
advected by gas motion, this preference in the direction of propagation is should be caused 
by supergranular flows. Figure [1] shows that the supergranular cell is somewhat elongated 
in X, suggesting that the flow is indeed stronger in that direction. 

We also compute the centrifugal acceleration vd(p/dt, where <fi is the direction of the 
IME motion. This is the relevant quantity when considering generation of transverse waves 
in fluxtubes. Measurements with velocities below 1 km/s are discarded, because the errors 
in the determination of the angle increases with decreasing velocity. The histogram of the 
centrifugal acceleration is shown in the right panel of Fig. [TOj The di stribution is Gaussian 



with an average of 0.045±0.013 and an rms of 4.559±0.011 deg km/s 2 . iNisenson et al.l (120031 ) 
found similar results, however, a detailed comparison is hampered by their low number of 
measurements. 
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Fig. 10. — Histograms of the direction of motion (left panel) and the centrifugal acceleration 
(right panel) of IMEs. Blue: positive polarity. Red: negative polarity. 




^ [deg] fraction of data 

Fig. 11. — Histogram of the distribution of angles between emerging bipoles (left panel), 
and the probability-probability plot of the distribution against the hypothesis of uniform 
distribution (right panel). The dashed line represents the cumulative distribution function 
of a uniform angular distribution. 

3.6. Emerging bipoles 

We search for emerging bipoles by pairing magnetic elements of opposite polarity that 
appear close in time and space. A pair of magnetic elements of opposite polarity is flagged 
if the elements appear within a radius of 1.4 Mm and no more than four time steps (140 s) 
apart. There are 639 pairs that satisfy these conditions. For these pairs we next compute 
the angle of the line connecting the center of the negative element with the center of the 
positive one. 

Figure [TT] shows the directional distribution and probability-probability ( "p-p" ) plot of 
the angle under the assumption that the angular distribution is uniform. A p-p plot can 
be used to see if a set of data follows a given distribution. It is constructed by plotting 
the cumulative distribution function F{xi) against (i — |)/n, where Xi is the ith data point, 
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ordered from smallest to largest, and i = 1,2, ... ,n. Here, we have taken a uniform distri- 
bution of angles F(xi) = (xi + 180)/360. If the data is perfectly uniformly distributed, we 
have Xi = 360 (i — |)/n — 180. Substitution yields F(xi) = (i — \)/n, so that a linear p-p 
plot results. Since the p-p plot in Fig. CD] is approximately linear, there does not appear to 
be preference in orientation of the emerging magnetic elements. We conclude that there is 
no evidence for Jo y's law at arcsecond scales from these data, in agreement with results of 



Lamb et all (120071 ). 



Pairs that are not actual emerging bipoles have random angles, and therefore only 
add a statistically uniform background that does not contribute signal to the directional 
distribution or the p-p plot in Fig. [TTJ However, the number of accidental associations may 
be large, so that any signal from real emerging bipoles is drowned i n noise. A sea r ch for 



emerging flux using Hinode's spectropolarimeter, such as the one by ICenteno et al.l (120071 ) 
would not suffer from potentially large numbers of false positives, even if automated, and 
would therefore give more trustworthy results. 



4. Summary and Conclusion 

We have analyzed the dynamics of IMEs using a sequence of magnetograms. Visual 
inspection of the data shows the existence of many long-lived magnetic elements that have 
frequent interactions with other elements during their lifetime. We find that they may mi- 
grate over distances of ~ 7 Mm in periods of several hours. Their interactions, migration, 
and the many short-lived concentrations of flux that appear in their vicinity make it cum- 
bersome to uniquely identify an IME, or a set of IMEs. The IMEs sometimes drop below 
the detection level in our data, but commonly reappear some time later. IMEs appear to 
outline cells on scales of several Mm. 

A manual inspection of IME locations in the G-band and Ca II H filtergrams shows that 
only about a fifth of the IMEs have associated bright points. Visual inspection reveals no 
obvious correlation between IME morphology and the existence of bright points. There is 
a substantial correlation between the existence of a bright point in the G-band and Ca II H 
intensities. Bright points are formed in the G band through weakening of molecular CH lines 
as a result of partial evacuation of the magnetic element, while they are caused in CallH by 
influx of radiation from the hot walls of the Wilson depression. We therefore attribute the 
similarity in the appearance of bright points in these passbands to their common origin. 

We identity magnetic elements using a hysteresis-based algorithm that is able to track 
IMEs for about 10 min on average. This is much shorter than the lifetime of several hours 
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predicted by iDe Wijn et aJj (120051 ). However, visual inspection shows that many elements 
intermittently drop below the detection limit of our algorithm, shortening their measured 
lifetime. 

IMEs exhibit proper motions that resemble a Gaussian distribution with a slight over- 
density of velocities near the origin and in the far wings. We measure an rms velocity of 
1.57 ±0.08 km/s. 

The IMEs show a slight bias of about 0.2 km/s for velocities toward the nearest network 
boundary. This bias persists up to timescales of at least 1200 s. It is the likely cause of weak 
positive velocity autocorrelation at long delay times. In the data analyzed here, IMEs of 
negative polarity show a statistically significant higher drift to the nearest network boundary. 
It is tempting to assume that this difference is somehow related to the dominant positive 
polarity of the network in the field of view. Alternatively, elements of negative polarity may 
be more mobile compared to elements of positive polarity because of their smaller size. 

We find a slight preference for velocities in X direction. Elongation in X of the main 
supergranular cell in the field of view suggests that supergranular flows are stronger in that 
direction. The observed preference in direction is thus in agreement with the classical picture 
that IMEs are advected by supergranular flows. 

IMEs experience centrifugal accelerations that obey a Gaussian distribution with an 
av erage of 0.045 ± 0.013 and an rms of 4.559 ± 0.011 deg km/s 2 , in agreement with results 



of iNisenson et all (12003 ). 



We search for emerging bipoles by pairing elements of opposite polarity that appear 
nearby each other in space and time. There is no detectable preference in the orientation of 
a pair. We conclude, therefore, that there is no evidence to support Joy's law at arcsecond 
scales from these data. 



Hinode is a Japanese mission developed and launched by ISAS/JAXA, collaborating 
with NAOJ as a domestic partner, NASA and STFC (UK) as international partners. Sci- 
entific operation of the Hinode mission is conducted by the Hinode science team organized 
at ISAS/JAXA. This team mainly consists of scientists from institutes in the partner coun- 
tries. Support for the post-launch operation is provided by JAXA and NAOJ (Japan), STFC 
(UK), NASA, ESA, and NSC (Norway). AdW thanks P. Judge, S. Mcintosh, R.J. Rutten, 
and J. Trujillo Bueno for discussions. 
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